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The problem of spectra formation in hydrodynamic approach to A+A 
collisions is considered within the Boltzmann equations. It is shown analyt- 
ically and illustrated by numerical calculations that the particle momentum 
spectra can be presented in the Cooper-Frye form despite freeze-out is not 
sharp and has the finite temporal width. The latter is equal to the inverse 
of the particle collision rate at points (i CT (r,p),r) of the maximal emission 
at a fixed momentum p. The set of these points forms the hypersurfaces 
t a (r,p) which strongly depend on the values of p and typically do not en- 
close completely the initially dense matter. This is an important difference 
from the standard Cooper-Frye prescription (CFp), with a common freeze- 
out hypersurface for all p, that affects significantly the predicted spectra. 
Also, the well known problem of CFp as for negative contributions to the 
spectra from non-space-like parts of the freeze-out hypersurface is naturally 
eliminated in this improved prescription. 

PACS numbers: 25.15.-q, 2^.10.Nz 

1. Introduction 

The Landau hydrodynamic approach [1] for multi-hadron production in 
hadronic/nuclear collisions appeared as a method that is alternative to the 
S-matrix one: the latter transforms the asymptotic hadronic states from 
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t = — oo to t = +oo while the former deals with the space-time evolution 
of thermal matter, produced in the collisions. The evolution is described 
by using the local energy-momentum conservation laws and thermodynamic 
equation of state of the matter. Unlike the S-matrix formalism, the initial 
state in hydrodynamic approach is associated with some concrete finite time 
after collision when created particles reach the locally equilibrated state and 
so, the initial state can be described by the minimal set of parameters. As 
well, because system formed is quite small, the picture of continuous medium 
is destroyed also at likely small finite time. Since the system expands fast, 
the latter is not constant in configuration space but depends on the position 
r of the fluid elements: t CT (r). The set of these points (i CT (r),r) in the 
Minkowski space forms, therefore, the hypersurface a corresponding to the 
outer boundary of the applicability of hydrodynamics. The sudden freeze- 
out implies that the spectra are formed just on this hyperurface where, 
formally, an ideal fluid transforms into an ideal gas. The particle momentum 
spectra then can be expressed by the well-known Cooper- Frye formula [2J: 



The freeze-out hypersurface is typically associated with an isotherm. Ac- 
cording to Landau PQ T ~ m n ~ 140 MeV. In current analysis of A+A 
collisions the corresponding temperature is in the region T = 90 — 150 MeV 
and is fixed, typically, from the best fit of the spectra. 

Any isotherm contains usually non-space-like parts, which lead to un- 
physical negative contributions to spectra for the particles with momenta 
directed inward the system, p^da^ < 0. There is no common phenomenolog- 
ical prescription, based on Heaviside step functions 9(p fM da /1 ), which allows 
one to eliminate the negative contributions to the momentum spectra when 
p^da^ < 0. The prescription, proposed in [3], eliminates the negative con- 
tributions in the way which preserves the number of particles in the fluid 
element crossing the freeze-out hypersurface. Therefore it takes into account 
that at the final stage the system is the only holder of emitted particles. An- 
other prescription [4] ignores the particle number conservation considering 
decaying hadronic system rather as a star - practically unlimited reservoir 
of emitted photons/particles. The both prescriptions have a problem with 
momentum-energy conservation laws at freeze-out. 

But the most serious problem is the obvious conflict of the CFp with 
simple observation: to provide sudden transformation of the liquid to ideal 
gas one needs to switch suddenly cross-sections between particles from very 
big values to very small ones [5] that cannot happen in reality. The rea- 
son of particle liberation is another: it is the gradual change of the ratio 
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between the rate of system expansion and the rate of collisions, so that 
particles should be emitted continuously. The phenomenological model of 
continuous emission was proposed in Ref. [6j; another approach is the so- 
called hybrid models [7] , where kinetic evolution is matched with hydro one 
at the hypersurface of hadronization. The shortcomings of these models 
are described in detail in [8j. But the feeling that the Landau/Cooper- Frye 
sudden freeze-out does not describe the real process of continuous particle 
emission is predominant now. 

Nevertheless, in Ref. [9] it was advanced an idea of duality in hydro- 
kinetic approach to A+A collisions: though the process of particle libera- 
tion, described by the emission function, is non-equilibrium and gradual, 
the observable spectra can yet be expressed by means of the Cooper-Frye 
prescription based on locally equilibrium distribution function. The con- 
clusive step towards an analytical and numerical realization of this idea is 
done in Ref. [8] where hydro-kinetic approach is developed: it overcomes all 
above mentioned problems by considering the continuous dynamical freeze- 
out that is consistent with Boltzmann equations and conservation laws. In 
what follows we stick to the mainstream of arguments developed in this 
paper. 



2. Kinetics of freeze-out in Boltzmann approach 

Let us start from Boltzmann equation. It has the general form: 

^^%A = Gi^ P )-L i ( X ,P). (2) 

The expressions G{(x,p) and Li(x,p) = R- L {x , p) fi(x , p) are so-called 
(G)ain and (L)oss terms for the particle of species i. Typically, i?j is a rate 
of collisions of z-th particle. Below we will omit index i, the corresponding 
expression can be related then, e.g., to pions. 

The probability Vt^>t'{ x iP) f° r a particle to reach the point x' = (f,r') 
starting from the point x = (t, r) without collisions is 

V t -+t'(x,P) = exp f- J dtR(x t ,p)^ , (3) 

where 

x t = (t,r + ^(t-t)). 

In terms of this probability, the Boltzmann equation can be rewritten in the 
following integral form 

f{t,r,p) = f(t ,r - -j-(t - t ),p)Vt ^t(t ,r fr(i - t ),p) 
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+ jG(T,r-^(t-T),p)V T ^t(r,r-^(t-T),p)dr. (4) 

to 

Let us integrate the distribution over the space variables to represent 
the particle momentum density at large enough time, t — > oo, when particles 
in the system stop to interact. To simplify notation let us introduce the 
escape probability for the particle with momentum p in the point x = (t, r) 
to leave system without collisions: V(t,r,p) = Vt-+( T -+Qo)(ti r iP)- Then the 
result can be presented in the general form found in Ref. [5]: 

n(t oo,p) = n(p) = J d 3 rf(t ,r,p)V(t ,r,p) 

oo 

+ J ' d 3 rJdt'G(t',r,p)V(t',r,p). (5) 
to 

The first term in Eq. ([5]) describes the contribution to the momen- 
tum spectrum from particles that are emitted from the very initial time, 
while the second one describes the continuous emission with emission den- 
sity S(x,p) = G(t,r,p)V(t,v,p) from 4D volume delimited by the initial 
and final (where particles stop to interact) 3D hyper surf aces. 

In what follows we will use the (generalized) relaxation time approx- 
imation proposed in [5], which is the basis of the hydro-kinetic approach, 
described in detail in [8]. Namely, it was argued [5] that there is such a local 
equilibrium distribution function fi_ eq _(T(x), u v {x), fJ-(x)) that, in the region 
of not very small densities where term G ~ S gives noticeable contribution 
to particle spectra, the function / is approximately equal to that one which 
would be obtained if all functions in r.h.s. of Eq. @ calculated by means of 
that function fi_ eq ,. The function fi eq ^ is determined from the local energy- 
momentum conservation laws based on the non-equilibrium function / in 
the way specified in [8]. Then, in accordance with this approach we use 

R(x,p) W R Leq Xx,p),G W Rl.eq.(x,P)fl.eq.(x,p). (6) 

The "relaxation time" r re i = 1/Ri. eq . grows with time in this method. 

3. Saddle point approximation for momentum spectra 

Let us generalize now the Landau/Cooper- Frye prescription (CFp) of 
sudden freeze-out. For this aim we apply the saddle point method to cal- 
culate the integral in the expression for spectra © with account of ©. To 
simplify notation we neglect the contribution to the spectra from hadrons 
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which are already free at the initial thermalization time to ~ 1 fm/c and 
thus omit the first term in 

To provide straightforward calculations leading to the Cooper-Frye form 
let us shift the spacial variables, r' = r+— (to— t 1 ), in ([5]) aiming to eliminate 
the variable t' in the argument of the function R which is the integrand in 
V(t',r,p). Then 

oo 

n(p) « J d 3 r' J dt'f^Xt'y + ?-(t' - t ),p)Q(t',v',p), (7) 
to 

where 

Q(t', r',p) = R(t', r' + ?-(t' - t ),p) exp { - 7 U( S) r' + ^-(s - *o),p)(fa 1(8) 

po ^ y po j 

Note that 

Q(t / ,r',p) = | 7 P(t',r',p), (9) 

where P(t',r',p) is connected with the escape probability V: 

P(t',r', P )=V(t',r' + ?-(t'-t ),p). (10) 
Po 

Therefore 

oo 

| dt'Q(t',r',p) = 1 - P(t ,r',p) « 1. (11) 

to 

The saddle point t a (r,p) is defined by the standard conditions: 

dQ{t',v>,p) 

-\t>=t> = U, 



dt' 

d 2 Q(t',r',p) 



-t> < 0. (12) 



dt* lt,: 

Then one can get from Q, (|10p the condition of the maximum of emission: 

P^dnR{t', T,p) ,11 P// 

~ R( t ' a ,r,p) \t'=t'„r=r'+^(t' a -t )=P0R(t (7 ,V + —(t a -t ),p). (13) 

If one neglects terms p*d r *R in l.h.s. and supposes that in the rest frame 
(marked be asterisk) of the fluid element with four-velocity u(x) the collision 
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rate, R*(x,p) = Po ^}^ ,p ^ , does not depend on particle momentum: R*{x) k, 
(v* a) (x)n* (x) (here n(x) is particle density, a is the particle cross-section, 
v is the relative velocity, < ... > means the average over all momenta), then 
the conditions (|13|) are equivalent to the requirement that at the temporal 
point of maximum of the emission function the rate of collisions is equal to 
the rate of system expansion [8]. This is the heuristic freeze-out criterion 
for sudden freeze-out [TO]. However, as we will demonstrate, the neglect of 
momentum dependence leads to quite significant errors. 

To pass to the Cooper-Frye representation we use the variables which 
include the saddle point: 

r = r'+ P-(t;(r»-t ). (14) 
Po 

Then the expression for the spectrum takes the form: 

oo 

n(p) « [d 3 r 1 - / dt'S(t',r,p), (15) 

J po or J 

to 

where the emission density in saddle point representation is [t a = (t a (r,p)) 

S(t',r,p) = /i.eq.(*',r + —(t' - t a ),p) 

Po 

xR{t a ,r,p)V(t a ,r,p) exp(-(i' - t a f/2D 2 (t a , r,p)). (16) 

According to Eq. ([3]) V(t a , r,p) = e , since the freeze-out zone is the region 
of the last collision for the particle. Then the normalization condition for 
Q (Q is presented by the bottom line in (|16p ) allows one to determine the 
temporal width of the emission at the point (t a (r,p), r,p): 

e 1 

D(t a ,r,p) = -=— « T iel (t a ,r,p). (17) 

V27T R{t a ,r,p) 



Therefore if the temporal homogeneity length A(i, r,p) of the distribu- 
tion function /i. eq . near the 4-point (i CT (r, p), r) is much larger than the 
width of the emission zone, X(t a , r,p) ^> r rel (t (T , r,p), then one can approx- 
imate /i.eq.(t',r + ^(t' - t a ),p) by /i. eq .(*a,r,p) in Eq. ([15]) and perform 
integration over t' accounting for normalizing condition (|11|) . As a result 
we get from (|15[) and (|16p the momentum spectrum in a form similar to the 
Cooper-Frye one ([!]): 



A(P)=/ W « J da^h. e<l Xx,p). (18) 

a(p) 
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It is worthy to note that the representation of the spectrum through 
emission function ([5]) is the result of the integration of the total non-equi- 
librium distribution function f(x,p), Eqs. ©, over the asymptotical 
hypersurface in time, while the approximate representation of the spectrum, 
Eq. (fT8|) . uses only the local equilibrium part /i. eq . of the total function 
f(x,p) at the set of points of maximal emission - at hypersurface (t a (r,p), r). 



Now let us summarize the conditions when the Landau/Cooper- Frye 
form for sudden freeze-out can be used. They are the following: 

i) For each momentum p, there is a region of r where the emission 
function as well as the function Q, Eq. ([8]), have a clear maximum. The 
temporal width of the emission D, defined by Eq. (|16p . which is found to 
be equal to the relaxation time (inverse of collision rate), should be smaller 
than the corresponding temporal homogeneity length of the distribution 
function: X(t a ,r,p) > D(t a ,r,p) ~ T iel (t a ,r,p). 

ii) The contribution to the spectrum from the residual region of r, where 
the saddle point method (Gaussian approximation (|16p and/or condition 
r rei "C A) is violated, does not affect essentially the particle momentum 
density. 

If these conditions are satisfied, then the momentum spectra can be pre- 
sented in the Cooper-Frye form despite the fact that actually it is not sudden 
freeze-out and the decoupling region has a finite temporal width T re i(t a ,r,p). 

The analytical results as for the temporal width of the spectra agree re- 
markably with the numerical calculations of pion emission function within 
hydro-kinetic model (HKM) [8j. For example, near the point of maximum, 
r = 16.5 fm/c, r = 0,pt = 0.2 GeV, the "experimental" temporal width 
Dhkm obtained by numerical solution of the complete hydro-kinetic equa- 
tions is Dhkm ~ 4.95 fm (see Fig. 1, left). Our theoretical estimate is 
D = J- 5.00 fm, since the rate of collisions in this phase-space point 

V 2ttH 



is R(T a (r,p) = 16.5 fm/c, r = 0,p T = 0.2 GeV, p L = 0) m 0.217 c/fm. 

It is worthy to emphasize that such a generalized Cooper-Frye represen- 
tation is related to freeze- out hyper surf aces that depend on the momentum 
p and typically do not enclose the initially dense matter. In Fig. 1, one 
can see the structure of the emission domains for different px in HKM |8j 
for initially (at r=l fm/c) Gaussian energy density profile with e max = 6 
GeV/fm 3 . The maximal emission regions for different pr are crossed by 
isotherms with different temperatures: 80 MeV for low momenta and 135 
MeV for high ones. This is completely reflected in the concave structure of 
the transverse momentum spectrum as one can see in Fig. 2. 



4. Generalized Cooper-Frye prescription 
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If a part of the hypersurface t a (r,p) is non-space-like and corresponds 
to the maximum of the emission of particles with momentum p, directed 
outward the system, the same part of the hypersurface cannot correspond 
to the maximal emission for particles with momentum directed inward the 
system. It is clear that the emission function at these points is close to zero 
for such particles. Even formally, in the Gaussian approximation (|16p for Q, 

oo 

validated in the region of its maximal value, the integral J dsR(s, r + ^(s — 

t a p 

t a (r,p),p) 3> 1, if particle world line crosses almost the whole system. The 
latter results in Q — > and, therefore, completely destroys the saddle-point 
approximation (|12p for Q and then the Cooper-Frye form (|18p for spectra. 
Recall that if a particle crosses some non-space-like part of the hypersurface 
a moving inward the system, this corresponds to the condition p^da^ < 
[3J. Hence the value p^da^p) in the generalized Cooper-Frye formula (|18|) 
should be always positive: p^da^(p) > across the hypersurface where fairly 
sharp maximum of emission of particles with momentum p is situated; and 
so requirement p^da^p) > is a necessary condition for t a (r,p) to be a 
true hypersurface of the maximal emission. It means that hypersurfaces of 
maximal emission for a given momentum p may be open in the space-time, 
not enclosing the high-density matter at the initial time to, and different 
for different p. All this is illustrated in Fig. 3, where the structure of 
particle emission domain is shown for two groups of particles. In the first 
one, the momentum is directed as the radius vector to the point of particle 
localization (they move outward the system), in the second one - in opposite 
direction (they move inward). The points of maximum for different pt, 
where Cooper-Frye form can be applied, do not overlap. The calculations 
have been done in HKM |8j . 

Therefore, there are no negative contributions to the particle momentum 
density from non-space-like sectors of the freeze-out hypersurface, that is a 
well known shortcoming of the Cooper-Frye prescription [3j[l]; the negative 
contributions could appear only as a result of utilization of improper freeze- 
out hypersurface that roughly ignores its momentum dependence and so is 
common for all p. If, anyhow, such a common hypersurface will be used, e.g. 
as the hypersurface of the maximal particle number emission (integrated 
over p), there is no possibility to justify the approximate expression for 
momentum spectra similar to Eq. (|18p . 

5. Conclusions 

Our analysis and numerical calculations show that the widely used phe- 
nomenological Landau/ Cooper-Frye prescription for calculation of pion (or 
other particle) spectrum is too rough if the freeze-out hypersurface is con- 
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sidered as common for all momenta of pions. The Cooper-Prye formula, 
however, could be applied in generalized form accounting for direct momen- 
tum dependence of the freeze-out hypersurface o~(p); the latter corresponds 
to the maximum of emission function S(t a (r,p),r,p) at fixed momentum 
p in an appropriate region of r. If such a hypersurface o~{p) is found, the 
condition of applicability of the Cooper- Frye formula for given p is that the 
width of the maximum, which in the simple cases - e.g., for one compo- 
nent system or at domination of elastic scatterings - is just the relaxation 
time (inverse of collision rate), should be smaller than the corresponding 
temporal homogeneity length of the distribution function. 
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Fig. 1. The pion emission function for different pt in hydro- kinetic model (HKM) 
[5]. The isotherms of 80 MeV (left) and 135 MeV (right) are superimposed. 
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Fig. 2. Transverse momentum spectrum of 7r in HKM, compared with the sudden 
freeze-out ones at temperatures of 80 and 160 MeV with arbitrary normalization. 
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Fig. 3. The emission function in HKM for particles with momentum directed along 
the radius vector at the emission points (left) and for those ones in the opposite 
direction to the radius vector (right). 



